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Abstract 

We explain the physical nature of the subset solution to the sign problem in chiral random 
matrix theory: The subset sum is shown to project out the canonical determinant with zero quark 
charge from a given configuration. As the grand canonical chiral random matrix partition function 
is independent of the chemical potential, the zero quark charge sector provides the full result. 
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I. INTRODUCTION 

Chiral random matrix theory |l| has given us several deep insights into the QCD sign 
problem whichprohibits direct application of lattice QCD methods at nonzero quark chem- 
ical potential [2[ . The random matrix framework has allowed us to understand the failure of 
the quenched approximation 3|, to formulate the OSV relation which replaces the Banks- 
Casher relation [5| and to derive the first analytical result for the average phase factor of 
the fermion determinant Q|. 

These lessons from chiral random matrix theory apply directly to QCD at nonzero chem- 
ical potential since the two are equivalent in the microscopic limit or the e-regime of chiral 
perturbation theory. For these reasons it is most interesting that the signproblem in a chiral 
random matrix theory can be solved by means of the subset method 7|, |8| . In particular 
this subset method works even in the region of fi > m n /2 where the sign problem is severe. 

The aim of the present paper is to provide the physical explanation of why the subset 
method introduced in [7, 8] solves the sign problem in chiral random matrix theory: As we 
will show in detail below the subset construction projects out the canonical determinant with 
zero quark charge from the fermion determinant. Since the chiral random matrix partition 
function is independent of the chemical potential the zero charge part makes up the full 
result. 

In this paper we start from a random matrix theory for QCD that is /i-independent even 
for finite size of the random matrix. This choice has a direct physical motivation: First, 
in the microscopic domain (where the size, n, of the random matrix goes to infinity while 
the quark mass times n and the square of the chemical potential times n are held fixed) the 
random matrix partition function is identical to the partition function of chiral perturbation 
theory in the e-regime 9] . Second, being a theory of pions, which are bound states of quarks 
and anti-quarks, chiral perturbation theory naturally does not couple to the quark chemical 
potential. As chiral perturbation theory is the effective theory for QCD at low temperatures 
and /i < m7v/3 the /z-indenpendence of the partition function is natural in this regime. 

In general, random matrix partition functions for QCD only need to be /i-independent 
in the microscopic domain. In the present context it is convenient to work with a random 
matrix theory where this /i-independence is manifest even at finite n. 
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II. SUBSET AND CANONICAL DETERMINANTS IN CHRMT 



Our starting point is a variation of the chiral Random Matrix Theory (chRMT) at non 



zero chemical potential [i introduced in 10| (see 



Z(m,fjL)= / rf$irf$ 2 det 



m e^i - e" M $ 
_ e -M$t _|_ e nQ 2 m 



11 1 for a review). It is defined by 



where $1 and $2 are complex n x n matrices. We have chosen to work with this form of the 
partition function because it is independent of \x even for finite n (this was also the case for 
the partition function used in 10(), and because the chemical potential appears in the form 
exp(±/i) which allows us to project out the canonical partition function in the same way as 
in lattice QCD. The /i-independence of the partition function follows immediately by using 
that the Gaussian integral is only nonzero for terms that have an equal number of factors 
$i and $| for % — 1, 2. The relation to the form used in 7, 8] is given in the appendix. 

In the subset method of {7, 8], one first performs a sum over a subset of roots of unity 
contained in the integral over the matrices $1 and $2- The critical observation in [7, 8| that 
the determinants 



/ 



d(fi, 9k) = det 



m 



_ e -M-* e fe$t _|_ e V+idk(§) 2 

where 9 k = 2kn/N s with N s > 2n + 1 



(2) 



m 



sum up to a positive real number 

- d (^ d k) e R+. 



N a -1 



N. 



(3) 



k=0 



This number can then in turn be used to generate a Monte Carlo ensemble of configurations 
and subsequently the unquenched expectation values. Note the invariance of the Gaus- 
sian measure under these phase rotations (the arguments below apply to any measure with 
the same invariance properties). Next we show that the measure Eq. (JSJ) is a canonical 
determinant with zero baryon number which is manifestly positive. 

For a given configuration ($1, $2) we decompose the fermion determinant 



D(fji) = det 



m 



e M $i — e" M $ 



(4) 



m 
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into canonical determinants 



2d 



where 



q=—2n 



(5) 



(6) 



(See 12l4l4J] for applications of canonical determinants to lattice QCD.) Likewise we decom- 
pose the partition function, Z, into canonical partition functions 

2n 
q=—2n 

where ((...) is the expectation value with respect to the Gaussian weight for $1,2) 



Z q = (D q ). 



(8) 



As Z is independent of /1 we necessarily have Z q = for q ^ 0. For odd q the canonical 
determinants vanish as well, D q= 2i+i = 0. This follows trivially from D(i(fi-\-7r)) = D(ifi) and 
exp(— iq(fi+ir)) = exp(— iqfi)(— l) q . For even index, however, the canonical determinants are 
nonzero for a typical configuration <3> 2 ), and only after averaging will one find Z q=2 i = 
for I ^ 0. 



To make the connection to the subset construction of 
partition functions 



3, |8| we first rewrite the canonical 



D q = — [ d9 e- iq{ - ifl+e) D(i(-ifi + 0)) 



(9) 

where in the first line we shifted the contour into the complex plane. The determinant inside 
the integrand is now 



D(fj, + i0) = det 



m 



(10) 



m 



To establish the relation between this subset construction and the canonical determinants 
introduced above first note that we can replace the subset-sum over Q k in Eq. by an 
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integral 

N' 



l N f^d^e k ) = Y f d9 dfao). (n) 

s fc=o 71 ^~ 7T 

This follows from the observation that the integrand is a polynomial in e ±ld of maximum 
order 2n and that all integrals follow from the orthogonality relations 

f f d9 e^ 6 = 5 jh (12) 

The same orthogonality relation holds for the sum 

AT S -1 

ir E ^ = h ( 13 ) 

S q=0 

provided that \j\, \l\ < N s . Therefore the sum over k gives the exact value of the integral if 
N s >2n + l. 

By comparison of Eq. (J2]) with Eq. f lTUj) we then see that the subset sum is equivalent to 
the projection onto the q = canonical determinant, that is 

, JV.-l 

— J2 dfa Ok) = D . (14) 

s k=0 

This is the physical explanation of what the subset is. 

Since Z q=0 = Z (the q = part makes up the entire partition function because it is 
independent of fi) the subset method gives the full result. Moreover, as the subset sum for 
a given configuration is equivalent to the canonical determinant with q = 0, it is clear that 
the subset sum is necessarily real and positive: as can be seen explicitly from Eq. ([9]), we 
have that D q= o is independent of fi, and for \i = all determinants in the subset sum are 
real and positive. This is the physical explanation of why the subset method works, see also 



For the variant of the chiral random matrix partition function used in |7|, |8| the interpre 
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tation of the subset is analogous, see Appendix A The original argument for why the subset 
method works given in [7, 8] is also related to the argument given above. 

In general the QCD partition function will of course depend on the chemical potential and 
hence in QCD one will need to evaluate all D q . For the evaluation of D q it is also possible 
to turn the integral into a sum. In this case the maximum order of the polynomial in e ±ld is 
2n + \q\, and therefore the subset sum evaluates the integral exactly for N s > 2n + \q\ + 1. 



Note, however, that the D q with q ^ are not real and positive so we do not have a weight 
to perform Monte Carlo Simulations (even if we can do these integrals exactly). This is in 
exact analogy with the observations of in lattice QCD. 

The argument given above also applies if $ is unitary rather than complex. More gen- 



erally, for unitary 
temporal links by 



attice gauge theories where the chemical potential is introduced into the 



15| 



Ut -> e^U\. (15) 

the partition function is /i-independent and equal to the charge zero canonical partition 
function. The subset method then applies in exactly the same way as in the random matrix 
model discussed above. 



III. CONCLUSIONS 

The subset solution to the sign problem in chiral random matrix theory has been shown to 
be equivalent to the projection, configuration by configuration, onto the zero quark number 
canonical determinants. Since the chiral random matrix partition function is independent 
of the chemical potential, the canonical partition function makes up the full grand canonical 
partition function. This gives the physical reason how the subset construction works. The 
same argument applies to unitary lattice gauge theories at nonzero chemical potential. 

The vanishing value of the canonical partition functions in chiral random matrix theory 
for nonzero quark number is the result of detailed cancellations: the canonical determinants 
with nonzero quark number take complex values and only the average value is zero. The 
projection onto the canonical determinants with nonzero quark number can also be obtained 
from a subset sum, however, it remains a challenge to devise a numerical method to control 
the cancellations in the average. Such a method would potentially have direct application 
to full QCD where partition functions with q ^ are nonvanishing. It may also be able to 
cast further light on the special nature of the noise [if], related to the sign problem. 

Despite the /i-independence of the chiral random matrix partition function the random 
matrix theory gives a plethora of nontrivial results for the spectral correlation functions 
of the Dirac operator and for the fluctuations of the fermion determinant. The reason for 
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this is that the generating functionals for such partially quenched observables have a highly 
nontrivial dependence on the chemical potential. It would be most interesting if one would 
be able to extend the subset method to these partially quenched observables. 
Acknowledgments: We would like to thank Christof Gattringer for discussions and ex- 
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Appendix A. EQUIVALENCE OF THE CHRMT FORMULATIONS 



The form of the chRMT used in 7J, [8j was 



'B{m,fi B ) 



d$id$ 2 det 



m 



z$i + // B $ s 



-nTr(*i3>t+$2 < I > 2) 



(16) 



m 



This partition function depends on \i B for finite n [10j. Because of the /independence of the 
partition function Z B (m, /is), the subset sum is not equal to the canonical partition function 
for q B = and the corresponding canonical partition functions for q B ^ are nonvanishing. 
The /^-dependence is, however, of a form where the partition function at non-zero fi B is 
trivially related to the one at /j, b — 10] 



Z B (m,fi B ) = (1 - n B ) n Z B ( 



m 



/4 



0). 



(17) 



In |8j it was shown that the subset sum for each configuration realizes this relation. When 



\i B < 1 both the prefactor 
thus, as originally argued in 



— [J> B ) n and the rescaled quark mass are real and positive 
7j, |8[ , the subset sum for the right hand side is always real and 
positive. The relation, Eq. (JT7J) , is the analogue of the \x- independence of the chRMT used 
in this paper, and the fact that subsets realizes this relation configuration by configuration 
is the analogue of the projection onto the canonical determinant with zero quark charge. 

The form of the chRMT adopted in Eq. (TfJ) is related to the form, Eq. ffl6|) . used in |7J, 
by a fi dependent rescaling of the mass and a trivial overall factor. If we start from 



Z(m,p) 



d$i<i<I>2 det 



my/l - jl 2 i$i + /i$ 2 
i${ + my/l - /t 2 



(18) 



then it is clear from Eq. (|T7|) that Z(m,jl) is independent of fi. Moreover, with \x given by 



tanh(/i) = fj,, (19) 
then this partition function is identical to the one of Eq. ([T]) . In order to see this first note 



that cosh(/i) = — /i 2 and sinh(/i) = /i/yl — /i 2 and then use this to express the 

determinant in terms of yU 



det 



m/ cosh(/i) 1/ cosh(/^)(i cosh(/i)$i + sinh(/i)$ 2 ) 

1/ cosh(/i)(z cosh(/x)<&J + sinh(/i)$|) mj cosh(/x) 



l/co S h-Mdet| m ico S h W * 1 + S i„h M $ 2 ^ (2Q) 

z cosh(/i)$J + sinh(/x)$2 m 

The factor 1/ cosh 2n (/i) cancels against the prefactor 1/(1 — jl 2 ) n in the partition function 
of Eq. (fl~8l) . After choosing 

^ = 1^-^), $' 2 = ^($l-z$t) (21) 

as new integration variables, we recover the form given in Eq. ([]]). 

nn 

The subsets defined in [3, [8| consist of rotated matrices 

$i ->■ cos^^i + sin ^$2, $ 2 - sin^^i + cos6 l A .$2 (22) 
which translates into 

$' 12 ^ e^$' 12) (23) 

as in Eq. (J2J. 
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